Accepted by ApJ 

Preprint typeset using I^T^jK style emulateapj v. 5/2/11 



TRACING GALAXIES THROUGH COSMIC TIME WITH NUMBER DENSITY SELECTION 

Joel Leja^ , Pieter van Dokkum^ , Marijn Franx^ 
Accepted by ApJ 

ABSTRACT 

A central challenge in observational studies of galaxy formation is how to associate progenitor 
galaxies with their descendants at lower redshifts. One promising approach is to link galaxies at 
fixed number density, rather than fixed luminosity or mass. This method is effective if stellar mass 
rank order is broadly conserved through cosmic time. In this paper, we use the Guo ct al. (2011) 
semi-analytical model to analyze under what circumstances this assumption is valid in the context 
of a cosmological simulation. Specifically, we select progenitor galaxies at a constant number density 
and compare the stellar mass evolution of their descendants to the evolution at a constant number 
density. The median stellar mass of the descendants increases by a factor of four (0.6 dex) from z = 3 
to 2: = 0. Constant number density selection reproduces this to within 40% (0.15 dex) over a wide 
range of number densities. We show that the discrepancy primarily results from scatter in the stellar 
mass growth rates and merging. After applying simple, observationally-based corrections for these 
processes, the discrepancy is reduced to 12% (0.05 dex). We conclude that number density selection 
can be used to predict the median descendant mass of high-redshift progenitor galaxies. The main 
uncertainty in this study is that semi-analytical models do not reproduce the observed mass evolution 
of galaxies, which makes the quantitative aggregate effects of star formation, merging, and quenching 
on the rank order of galaxies somewhat uncertain. 

Subject headings: cosmology: dark matter — cosmology: theory — galaxies: abundances — galaxies: 
evolution — galaxies: statistics 



1. INTRODUCTION 

In order to infer changes in a galaxy population, it is 

necessary to link progenitor and descendant galaxies at 
different snapshots in redshift. Yet in order to properly 
link these galaxies, it is also necessary to know how their 
observational properties change through time i.e., it is 
necessary to already have a consistent theory of galaxy 
evolution. This creates a pseudo "chicken and egg" prob- 
lem which has been a long-standing challenge for studies 
of galaxy evolution. 

One common solution to the progenitor problem is to 
link galaxies at fixed stellar mass at different redshifts. 
This is effective when studying galaxies which are be- 
lieved to evolve passively, as their properties are rela- 
tively constant with time: for example, luminous red 
galaxies (LRGs) (Wake et al. 2006). This method, how- 
ever, suffers from progenitor bias (van Dokkum & Franx 
1996): by only studying galaxies at high redshift that 
appear similar to local LRGs, some of the progenitors 
of local LRGs will be missed. Linking galaxies in this 
fashion will thus systematically underestimate their true 
evolution. As the lookback time increases, the error in- 
trodiiced by comparing galaxies at fixed mass will also 
increase, making the method unsuitable for studies be- 
yond the nearby Universe. 

Some recent studies (Wuyts et al. 2011; Leitner 2012) 
suggest linking local starforming galaxies to their pro- 
genitors at high redshift by following the mass evolution 
implied by the starforming sequence. The tight relation- 
ship between stellar mass and star formation rate makes 
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this a powerful technique. However, for this techniqiie 
to work, the mass added via mergers must be negligible. 
Additionally, systematic uncertainties in the star forma- 
tion rate (SFR) are considerable and vary substantially 
between SFR indicators (Wuyts et al. 2011). The tech- 
nique also cannot be applied to quiescent galaxies, as 
their quenching redshift and subseqiicnt merger history 
remains unknown. Since many galaxies are quiescent 
(Brammer et al. 2011), other methods to link descen- 
dants and progenitors must be developed. 

Tracking galaxies at constant number density is a 
promising technique that may prove more effective than 
previous methods. The basic approach, introduced in 
van Dokkum et al. (2010), is to assume that the number 
density of galaxies when ordered by some physical prop- 
erty (e.g., stellar mass) does not evolve. This provides 
a simple way to link descendant and progenitor popula- 
tions at different redshifts. This technique has already 
been used to study the evolution of a number of differ- 
ent galaxy properties: UV luminosity and star formation 
history (Papovich et al. 2011), the stellar velocity dis- 
persion function (Bezanson et al. 2011), Ha equivalent 
width (Fumagalli et al. 2012), and mass and structural 
evolution (Bezanson et al. 2009; van Dokkum et al. 2010; 
Brammer et al. 2011; Patel et al. 2012). 

However, investigations into the effectiveness of the 
technique have been limited. The investigations that do 
exist rely primarily on numerical results. Van Dokkum 
(2010) used a simple Monte Carlo simulation in their 
appendix to show that mergers have little effect on the 
stellar mass evolution measured at a constant number 
density. Papovich et al. (2011) went slightly further, 
using halo merger trees from the Millenium Simulation 
(Springel et al. 2005) to show the recovery fraction of 
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true descendant halos from z~7 to Z'-^S is approximately 
50%. Dark matter halo mass evolves in a fundamentally 
different way than stellar mass, though, and can be sub- 
ject to substantial tidal stripping when a halo becomes 
a satellite: this makes it less suitable for galaxy tracking 
than stellar mass. No study to date has attempted to ex- 
amine the efficacy of number density selection in a more 
realistic simulation setting. Additionally, the theoretical 
underpinnings of number density selection remain unex- 
plored: the rules governing how it works, why it works, 
and where it breaks down remain an open question. 

This paper will study number density selection by first 
discussing the assumptions inherent in the use of number 
density selection and which physical processes might vi- 
olate them. The Guo et al. (20f 1) semi-analytical model 
(SAM) applied to the Millenium Simulation is then used 
to investigate the validity of these assumptions. Semi- 
analytical models represent our best theoretical formula- 
tion of how galaxy evolution works on a broad scale, and 
simulations are currently the only way in which it is pos- 
sible to track the actual evolution of individual galaxies. 
In addition, it is not necessary for our purposes that the 
SAM accurately reproduces the observed Universe. This 
makes the SAM useful as a simple tool to explore the me- 
chanics of number density selection and what processes 
may cause it to err. 

In Section 2, we introduce the number density selection 
technique and discuss potential sources of error. In Sec- 
tion 3, we describe the Millenium Simulation, the Guo 
et al. (2011) semi-analytical model, and describe how 
number density selection is applied to these data. In 
Section 4, we present the resulting median stellar mass 
evolution and several techniques to improve it. In Sec- 
tion 5, we explore the possibility of ordering galaxies by 
inferred velocity dispersion rather than stellar mass. In 
Section 6, the suitability of applying SAM-based results 
to the real Universe is discussed. The conclusions are 
found in Section 7. 

When necessary, we assume Ho=70 km s^^ Mpc"'^ 
throughout the paper. 



2. CONSTANT NUMBER DENSITY SELECTION 

We discuss the underlying assumptions and potential 
pitfalls of the number density selection technique here. 
As described earlier, linking progenitors and descendants 
at a constant number density implicitly assumes that the 
number density of galaxies does not evolve. One can 
then infer the mass evolution of a given galaxy popu- 
lation from the redshift evolution of the galaxy stellar 
mass function. Number density selection thus provides 
a physically-motivated, redshift dependent stellar mass 
selection. The preferred approach is to link galaxies at a 
constant cumulative number density, which has the ad- 
vantage of being single-valued in mass (Brammer et al. 
2011). Additionally, the physical interpretation is sim- 
ple: if the stellar mass rank order of galaxies remains the 
same through cosmic time, then tracking galaxies at a 
cumulative number density will be effective. 

The aim of this study, then, is to investigate the condi- 
tions under which the conservation of rank order holds. 
If stellar mass growth is either constant or monotonically 
increasing in proportion to stellar mass itself, then rank 
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Fig. 1. — Stellar mass assembly histories are shown for four repre- 
sentative galaxies from the Guo et al. (2011) semi-analytical model 
output. Despite similar initial stellar masses, their mass growth 
histories vary substantially. The primary mechanisms for stellar 
mass growth (major and minor mergers, star formation) are la- 
belled. Furthermore, the processes interfering with number density 
selection- scatter in growth rate, merged galaxies, and quenching- 
are also labelled. In order to for constant number density selection 
to produce robust results, these processes must be either negligible 
or correctable. 

order will be preserved. The stellar mas^of a galaxy can 
be increased in two ways: star formation and merging. 
Steady star formation in starforming galaxies is propor- 
tional to M" (where a > 0) over a wide range of red- 
shifts, commonly referred to as the starforming sequence 
(Brinchmann et al. 2004; Noeske et al. 2007; Peng et al. 
2010; Gonzalez et al. 2012; Whitaker et al. 2012). The 
galaxy-galaxy merger rate is also expected, on average, 
to increase monotonically with stellar mass (see, e.g., the 
review by Hopkins et al. (2010)). Thus, there is good 
reason to believe rank order in stellar mass may be con- 
served. 

However, there also exist physical processes which will 
confuse rank order. Measured scatter in the starforming 
sequence is 0.34 dex; this includes intrinsic scatter from 
AGN, starbursts, and galaxies in the process of shut- 
ting down their star formation as well as measurement 
error (Whitaker et al. 2012). The stochastic nature of 
mergers naturally induces scatter in stellar mass growth 
over short timescales. These two effects in combination 
are referred to as the scatter in the stellar mass growth 
rate, and can break rank order by either "jumping" a 
galaxy ahead or "leaving" a galaxy behind. An exam- 
ple of breaking rank order by scatter in growth rates is 
shown in Figure [l] by the green galaxy trajectory. 

Galaxy mergers will also change the cumulative num- 
ber density by decreasing the total number of galaxies 
above a fixed stellar mass. This will cause a constant 
number density line to underestimate the true median 
mass evolution. This issue is only relevant for merg- 

^ Note that stellar mass in the Guo et al (2011) model refers to 
the integral of the star formation rate. Guo et al. (2011) do not 
explicitly include mass loss from SNe and winds in their model; 
instead a mass fraction of 0.43 is immediately returned to the ISM. 
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ers with galaxies that are at least as massive as the 
galaxy under consideration. As lower-mass galaxies have 
a higher number density than the number density under 
consideration, they do not affect the inferred stellar mass 
evolution. 

For completeness, one must also consider the creation 
of new galaxies. The stellar mass at which this effect is 
important is well below current detection limits in high- 
redshift surveys. For example, the mass-complete limit 
for the Newfirm Medium Band Survey (Whitaker et al. 
2011) is M* > 10" Mq (all galaxies) or M* > lO^" Mq 
(starforming galaxies) at z=2.2 (Brammer et al. 2011). 
Assuming newly formed galaxies at this epoch lie on the 
starforming sequence, they would not approach the mass 
limit in any reasonable amount of time. The creation of 
new galaxies may then be safely neglected. 

Finally, quenching may significantly affect rank order. 
Quenched galaxies no longer grow via star formation, 
whereas their star-forming companions at a similar mass 
will continue to add stellar mass. For central galaxies, 
this may be offset somewhat as they still exp erie nce sig- 
nificant growth through mergers: see Section [4. 3| for fur- 
ther discussion. Redshift evolution of the quenching rate 
and quenching mass limit may also confuse the rank or- 
der of galaxies. 

The effects of these processes on the rank order of 
galaxies are illustrated in Figure [T] and must be well 
characterized in order to safely use number density se- 
lection to track galaxy evolution. 

3. SIMULATION DATA AND SELECTION 
3.1. The Millenium Simulation 

We test the number density selection technique using 
mock galaxy catalogs from the Guo et al. (2011) semi- 
analytical model, which is based on N-body output from 
the Millenium Simulation (Springel et al. 2005). We 
access these data using the online relational database 
built by the German Astrophysical Virtual Observatory 
(Lemson & Virgo Consortium 2006). 

The Millenium Simulation is a large cosmological N- 
body simulation following the evolution of 2160"^ parti- 
cles from z = 127 to the present. The Millenium cosmol- 
ogy is based on a joint analysis of the 2dFGRS (CoUess 
et al. 2001) and WMAPl results (Spergel et al. 2003). 
The input cosmological parameters are JIa = 0.75, f2„i 
= 0.25, Qb = 0.045, 71 = 1, erg = 0.9, and Hq = 73 
km s~^ Mpc~^. While most of these parameters have 
changed in the recent WMAP7 results (Komatsu 2010), 
the deviation is small enough to have little effect on the 
results (Guo et al. 2011). The only significant devia- 
tion is in CTg, which is almost 4cr off the WMAP7 value 
((78 = 0.809 ±0.024). The effect of an inflated erg parame- 
ter on mer ger rates and galaxy catalogs will be discussed 
in Section [6Tl 

The Millenium Simulation is run in a periodic box 
500h^^ Mpc on a side with a particle mass of 8.6x 10*h~^ 
Mq. Post-processing is performed to separate particles 
into fricnds-of- friends groups (Davis et al. 1985). The 
SUBFIND algorithm (Springel et al. 2001) is then used 
to split friends-of-friends groups into bound subhalos. 
The largest subhalo is referred to as the main subhalo, 
and contains most of the mass. A subhalo must contain 
at least twenty particles to be identified by the SUB- 




FlG. 2. — The stellar mass function for four redshift slices from 
the Guo et al. (2011) mock galaxy catalogs is illustrated here. The 
three number density slices investigated in this paper are shown as 
dotted lines. The intersection of the constant number density lines 
with the mass functions represent the expected mass evolution of 
the number density selected sample. This expected median mass 
evolution is also shown here as a function of number density. 



FIND algorithm, which corresponds to a minimum sub- 
halo mass of 1.72 x 10^° Mq. 

3.2. The Guo et al. (2011) SAM 

The Guo et al. (2011) SAM is an updated version 
of the earlier Milleniuni-based galaxy formation models 
(De Lucia et al. 2006). It successfully matches the abun- 
dance of '^lO^^ Mq galaxies out to z ~ 3. Additionally, 
it accurately reproduces the local stellar mass function 
by construction. However, due to the inflated trg param- 
eter in the Millenium cosmology, galaxies cluster more 
strongly than in the real Universe. This results in a high 
fraction of quenched red satellites. In order to match the 
local mass function, the SAM must assemble these low- 
mass quenched satellites early in the simulation. Thus, 
the abundance of low-mass galaxies beyond z = 0.6 is 
overpredicted. The effect of this on constant number 
density selection is discussed in Section [6j 

Additionally, due to the finite resolution of the N-body 
simulation, the mock galaxy catalogs are only reliable 
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TABLE 1 
Number, Density Samples 



n [Mpc^ 




An 




log(M*//i) [z=3] 


5 X 10- 


-5 


5 X 10- 


5 


10.48 


2 X 10" 


-4 


1 X 10- 


4 


10.31 


8 X 10" 


-4 


1 X 10- 


4 


10.02 



above a certain stellar mass limit. Application of the 
SAM to the higher-resolution Millenium II indicates that 
the mass functions and morphological classifications be- 
gin to diverge substantially around log Alstar ~ 9.8 at 
z = (Guo et al. 2011). Number density slices are cho- 
sen carefully so as to track galaxies only when above this 
stellar mass limit (see Section 3.3). Galaxy evolution 



subsequent to crossing the resolution limit is expected to 
behave with more regularity. 

3.3. Number Density Selection in the Guo et al. (2011) 

SAM 

Here we discuss the implementation of number density 
selection within the SAM. 

First, the mock galaxy catalog is used to build stellar 
mass functions at the redshift outputs nearest to z = 3, 
z = 2, z = 1, and z = (see Figure |2]). This covers 
the typical range of observational studies and includes 
the last '-^ll Gyr of galaxy formation. Galaxies at three 
different number densities are selected at z = 3, and 
their descendants are followed until z = 0. The mass 
evolution of these descendants is compared to the mass 
evolution inferred at a constant number density. We note 
that measuring the stellar mass at a fixed number density 
makes no use of a bin. 

Merging not only grows the stellar mass of galaxies, 
but also decreases the total number of galaxies. In order 
to self-consistently count galaxies, galaxies that merge 
with a more massive galaxy are considered to have "dis- 
appeared", and are discarded from the sample. As a 
consequence of this, the number of galaxies descended 
from the z = 3 number density bin will decrease as time 
pass es. T he magnitude of this effect is discussed in Sec- 
tion g^l 

The chosen number density slices are shown in Table 
1, along with the initial bin width at z = 3 and the corre- 
sponding mass at constant number density. The number 
density slices were chosen to sample a variety of stellar 
masses while remaining above the stellar mass resolution 
limit of XogMstar ~ 9.8 (see the discussion in Section 
3.2|). The bin size is constant in number density rather 



than in stellar mass for two reasons: (1) the median stel- 
lar mass is invariant with respect to bin size for a bin 
in number density, and (2) a bin size in number density 
provides a self-consistent way to compare bin size if one 
orders galaxies by a different parameter, as is done in 
Section [5] . Additionally, a bin size specified in units of 
cumulative number density means the projected number 
density of selected galaxies in constant in each redshift 
interval, which may be useful in observational studies. 

In order to ensure the results of this study are inde- 
pendent of the method of investigation, the bin size is 
varied for the n = 2 x 10~^ Mpc~'^ slice. The resulting 
effect on descendant mass evolution is shown in Figure 
[3] The descendant mass evolution is very stable with re- 
spect to the bin size. Also shown is the recovery fraction, 
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Fig. 3. — Examining the effect of varying bin size with the 
n = 2 X lO"* Mpc~'^ slice in number density. The recovery frac- 
tion is defined as the number of descendant galaxies within the 
bin divided by the total number of surviving descendant galaxies, 
while the contamination fraction is defined as the number of galax- 
ies in the bin not descended from the original selection divided by 
the total number of galaxies within the number density bin. Both 
are shown in the upper panel. The lower panel shows the change 
in the measured median stellar mass of the descendants. Bin size 
has very little effect on the recovered stellar mass, whereas recov- 
ery fraction and contamination fraction are strong functions of the 
chosen bin size. This suggests that recovery fraction and contam- 
ination fraction poor diagnostics to test number density selection, 
whereas median stellar mass is a much more robust indicator. 



defined as the number of descendant galaxies within the 
bin divided by the total number of surviving descendant 
galaxies, and the contamination fraction, defined as the 
number of galaxies in the bin not descended from the 
original selection divided by the total number of galax- 
ies within the number density bin. The recovery and 
contamination fractions are a function of (1) the median 
mass growth rate, (2) the dispersion in mass growth rate, 
and (3) the initial bin size. While these diagnostics have 
figured prominently in previous studies of number den- 
sity selection [e.g. van Dokkum et al. (2010); Papovich 
et al. (2011)], they are not discussed in this paper as 
they are demonstrated to be a strong function of bin 
size. Additionally, it is not necessary to recover the ex- 
act same galaxies with number density selection as long 
as the recovered galaxies are indistinguishable from the 
true descendants in important galaxy properties such as 
stellar mass. 

4. TRACKING STELLAR MASS EVOLUTION 

We follow the median stellar mass of galaxies from 
the original z = 3 sample and compare it to the stellar 
mass evolution at a constant cumulative number density. 
Galaxies that merge with a more massive galaxy are re- 
moved from the sample, as described in Section [XS] The 
results are shown in Figure [4] 

The median stellar mass of the descendants is repro- 
duced to within 0.15 dex over 0.6 dex of evolution in mass 
and ~ 11 Gyr in time for all number density slices. This 
corresponds to a factor of four growth in stellar mass 
predicted to within 40%. Excluding the z = point 
for the lowest number density bin, the error decreases to 
20%, or 0.08 dex. As time passes, the discrepancy be- 
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Fig. 4. — The median mass growth of the true descendants (sohd 
hne) compared to the mass growth inferred at a constant number 
density (dashed hne) for three bins in number density. The mass 
evolution is reproduced correctly to within 40% (0.15 dex) over an 
evolution of ~4 in mass (0.6 dex) and 11 Gyr in time. Comparing 
galaxies at the same stellar mass would be off by a factor of four 
(0.6 dex). Note that the mass evolution of the highest number 
density sample is systematically underpredicted, while the mass 
evolution of the two lower number density samples are systemati- 
cally overpredicted. This can be explained by the scatter in growth 
rates and the merging of galaxies, as described in Sections 4.1 and 
4.2, respectively. 

tween the predicted evolution and the actual evolution 
increases. This is because processes which scramble rank 
order have more time to operate. 

The direction of the discrepancy between the median 
descendant stellar mass and the stellar mass at a constant 
number density is consistent within each number density 
sample. The two lower number density samples over- 
predict the mass evolution of the descendants while the 
highest number density sample underpredicts the mass 
evolution of the descendants. The physical processes re- 
sponsible for this are identified in Figures [5] and [6j and 
corrected for in Figure [7] Discussion of these procedues 
is found in the following sections. 

4.1. Scatter in Growth Rates 

The two lower number density bins consistently over- 
predict the mass evolution of the sample. One plausible 
explanation for this is the scatter in stellar mass growth 
rates. As the stellar mass function is steep for these low 
number densities, there is a large reservoir of slightly 
lower mass galaxies. Due to the scatter in growth rates, 
some of these low mass galaxies will grow faster than 
the descendants of the z = 3 sample, overtaking them in 
stellar mass: this is a version of the well-known Edding- 
ton bias (Eddington 1913). To investigate this possibil- 
ity, the distribution of actual growth rates for individual 
galaxies in each number density bin is shown in Figure 

El . 

There is a peak at null growth in each number density 
bin. This peak is populated with galaxies that experi- 
enced no star formation nor significant growth via merg- 
ers over the time period. This peak is initially largest 
for the lowest number density population. However, as 
gas supplies dwindle, more and more galaxies become 
quenched. By z = 0, the highest number density bin 
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Fig. 5. — Tracing the evolution of descendants from the initial 
selection at z = 3.06. The panels show the distribution of evolution 
in stellar mass for individual galaxies. The average is marked with 
a circle, while the median is marked with a cross. The difference 
between the 75th quartile and the 25th quartile is shown in the 
upper left of each panel. The n = 5 X 10~^ Mpc~'^ descendants 
are multiplied by a factor of two for comparison, to compensate for 
the smaller initial bin size. 



has the largest number of galaxies that exhibit little to 
no growth. This is likely because the lower number den- 
sity bins contain more massive centrals which are still 
experiencing growth via merging. 

There is also a secondary peak in each growth rate 
which is most evident at the highest redshift. This is the 
peak due to steady star formation from galaxies on the 
star-forming sequence. Its width is a combination of the 
finite size of the mass bin and scatter in the star forming 
sequence. As more galaxies quench, this peak washes 
out. 

The median mass growth is marked with a cross. Ini- 
tially, the highest number density galaxies experience the 
largest relative growth. Note that this is not in violation 
of the assumption that more massive galaxies grow more 
quickly: in linear units, the lowest number density bin 
adds more mass than the higher number density bins. 
By the final redshift bin, as merging begins to dominate 
the stellar mass growth, this has reversed, and the high- 
est number density galaxies grow at a faster relative rate 
than the lower number density galaxies. 

An estimate of the effect of scatter in growth rates on 
number density selection can be made by convolving the 
stellar mass functions shown in Figure[2]with the individ- 
ual growth rates shown in Figure [5] At each snapshot in 
redshift, the growth rate for each number density slice is 
convolved with the stellar mass function at that redshift. 
This will "grow" each galaxy in mass with a probabil- 
ity distribution given by Figure [5j Due to the shape of 
the mass function, the convolved stellar mass function 
results in a larger inferred mass growth than the uncon- 
volved stellar mass function for all three number densi- 
ties. The difference between the mass growth inferred 
from the convolved mass function and that inferred from 
the original mass function is taken to be the error in in- 
ferred growth rate for each number density slice. The 
inferred mass growth adjusted for this effect is shown in 
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Fig. 6. — The cumulative percentage of galaxies that have lost 
their identity by merging into a more massive galaxy is shown. The 
percentage is plotted as a function of number density and redshift. 
The fraction of merged galaxies increases with number density, as 
galaxies at a lower number density have few galaxies more massive 
with which to merge. 



the bottom left panel of Figure [7] 

This procedure is only an approximation, as it assumes 
the measured growth rate within a number density slice 
is applicable to all galaxies in the SAM. However, this 
is a relatively reasonable approximation, as the galaxies 
that do "overtake" more massive galaxies are likely to be 
fairly close in stellar mass. As shown in Figure [Tj the 
correction successfully removes much of the discrepancy 
between the inferred stellar mass and actual stellar mass 
at the lower number densities. 

4.2. Merging Galaxies 

The highest number density bin, in contrast, systemat- 
ically underpredicts the stellar mass evolution of the true 
descendants. This effect becomes more pronounced after 
adjusting for the scatter in growth rates. A likely expla- 
nation is the change in cumulative number density caused 
by galaxies that " disap pear" due to merger events, as de- 
scribed in Section |3.3| When a galaxy with stellar mass 
greater than the sample disappears in a merger, this will 
decrease the actual cumulative number density of the 
sample. This results in systematic underprediction of 
the stellar mass evolution. 

The magnitude of this effect can be calculated exactly 
if the rate of galaxy-galaxy mergers as a function of stel- 
lar mass is known. A simple estimate can be made by 
calculating the fraction of galaxies that disappear due to 
mergers as a function of redshift, shown in Figure [6j The 
cumulative number density is multiplied by the merger 
fraction. As the cumulative number density is decreased, 
the inferred mass growth will increase. The inferred mass 
growth corrected for merged galaxies is shown in the top 
right panel of Figure [7] 

This correction is only approximate, as it implicitly 
assumes that the merger fraction measured at a given 
number density holds for all galaxies at lower number 
densities. A more precise way to correct for this effect 
would be to integrate the merger fraction rate as a func- 
tion of number density. This would be a small adjust- 



ment for the higher two number density bins, as can be 
seen in Figure [6j since the merger fraction is relatively 
constant over number density in this regime. The lowest 
number density bin includes the most massive galaxies in 
the simulation, however, where the merger fraction must 
necessarily go to zero at the massive end. It is likely 
then that this approximation overadjusts the mass evo- 
lution in the lowest number density bin. The fact that 
the fully "corrected" inferred mass growth is larger than 
the actual mass growth lends support to this hypothesis. 

The combined effect of both corrections is shown in 
the lower right panel of Figure [7] The agreement is now 
excellent, with a discrepancy ofonly 12% (0.05 dex) over 
a factor of four (0.6 dex) evolution in stellar mass. 

4.3. Quenching 

In order for number density selection to be effective, 
the stellar mass growth of galaxies must be a monoton- 
ically increasing function of stellar mass. This is true 
separately for both modes of growth: star formation and 
merging. However, if star formation is shut down in a cer- 
tain subset of the population, a discrepancy may arise be- 
tween the mass growth rate of starforming and quenched 
galaxies. This will confuse rank order. We investigate 
the effect by using a color cut to separate quiescent and 
starforming galaxies at each redshift (see Figure [8| and 
then plotting the distribution of mass for both popula- 
tions as a function of redshift for the n = 2 x 10"'' Mpc~'^ 
number density slice (see Figure^. 



It is evident from Figure 
tematic difference between t 



'9^ that there is a clear sys- 
re mass distribution of star- 



forming and quiescent galaxies at most redshifts. The 
median mass of starforming descendants is higher than 
the median mass of quiescent descendants, except (by 
construction) the initial selection at z = 3. After the 
time between z = 3 and z = 2, however, the median mass 
growth of quiescent galaxies mirrors that of starforming 
galaxies. This offset is small (0.15 dex) compared to the 
spread in stellar mass of the full sample of descendants 
(full-width half-max 0.5 dex). However, it motivates 
a closer examination of the processes at work. 

A number of low-mass quiescent galaxies are estab- 
lished early on. Due to their low mass, they are likely to 
be satellites quenched by environmental processes. This 
low-mass tail is present throughout the simulation, al- 
though it is not necessarily composed of the same galax- 
ies in each snapshot. 

The number of quiescent galaxies grows by the quench- 
ing of star formation in starforming galaxies. Since star- 
forming galaxies in the sample have a higher median 
mass than quiescent galaxies, quenching will necessar- 
ily grow the upper envelope of the quiescent mass dis- 
tribution. However, at the highest masses, there exist 
quiescent galaxies that are more massive than any star- 
forming galaxies. Despite being quiescent, these galaxies 
must have substantially increased their stellar mass after 
quenching: the only logical avenue for this is growth via 
mergers. 

This leads to the possibility of three unique galaxy pop- 
ulations, distinguished by their method of stellar mass 
growth (or lack thereof): starforming galaxies, quenched 
centrals, and quenched satellites. This is illustrated fur- 
ther in Figure [To] The left panel plots star formation rate 
versus stellar mass, with populations split into quiescent 
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Fig. 7. — The median stellar mass of the true descendants and the corresponding stellar mass inferred from the stellar mass function is 
shown in three different number density bins. The top left panel shows the median stellar mass with no adjustments. The top right panel 
shows the median stellar mass after correcting the cumulative number density for galaxies that have been destroyed via mergers. The 
lower left panel shows the median stellar mass after adjusting the number density selection to accommodate for the scatter in stellar mass 
growth. The final panel shows the combination of these two corrections. The remaining error in the inferred mass evolution is likely due 
to a combination of the effects of quenched satellites and the approximate nature of the corrections. 



I 0.4 




9.5 10.0 10.5 11.0 11.5 
stellar moss [log(MQ/h)] 



12.0 



Fig. 8. — The color-mass relation at z = in the Guo et al. 
(2011) semi-analytical model is shown here as an example. The 
division between starforming and quiescent galaxies is indicated by 
the blue dashed line. At higher redshlft, this division is adjusted 
accordingly. 



(red) and starforming (blue) by a color-mass cut. The 
running median is shown for both. Clearly, this presents 
an issue for number density selection, as the quiescent 
galaxies seem to be adding mass at a much slower rate. 
However, this only accounts for stellar mass added by 
star formation. 

The right panel of Figure [TO] includes stellar mass 
added by mergers (averaged over z = — 0.28) as well 
as star formation. This is the total mass growth rate as 



a function of stellar mass. The plot looks markedly dif- 
ferent; at masses above ~ 10^^ M©, the median mass 
growth rate for quiescent and starforming galaxies is 
approximately the same. It is plausible, then, that 
quenched centrals and starforming galaxies at the same 
stellar mass are adding stellar mass at the same rate, 
thus preserving rank order for these galaxies. Quenched 
satellites remain an issue, as they have ceased adding 
mass altogether. We caution that these results are sen- 
sitive to the details of the quenching model in the SAM 
and the cosmology of the Millenium Simulation. 

5. ORDERING GALAXIES BY INFERRED VELOCITY 
DISPERSION 

Throughout this paper we have used stellar mass to 
order galaxies. Another promising approach to number 
density selection involves ordering galaxies by their stel- 
lar velocity dispersion. Velocity dispersion is a measure- 
ment of orbital energy and gravitational potential depth. 
There is evidence from both observations and simulations 
that velocity dispersion may be a more stable quantity 
than stellar mass (Loeb & Peebles 2003; Franx et al. 
2008; van Dokkum et al. 2010; Bezanson et al. 2011; 
Oser et al. 2012). Additionally, velocity dispersion has 
been proven to be a better predictor than stellar mass 
of specific star formation rate and color (Franx et al. 
2008; Wake, van Dokkum, & Franx 2012b), galaxy clus- 
tering properties (Wake, Franx, & van Dokkum (2012a)- 
but see also Li, Wang, & Jing (2012)), and star forma- 
tion rate (Bezanson, van Dokkum, & Franx 2012) over 
a wide range of redshift. This motivates an attempt to 
track galaxies by their stellar velocity dispersion, rather 
than their stellar mass. 
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Fig. 9. — The left panel shows the mass distribution of the descendants of the n = 2 X 10"'' Mpc"'' number density panel, split into 
starforming and quiescent galaxies based on a color-mass cut (see Fig. 7). The y-axis of the inset plots is dN/dlog(A/siar) with a consistent 
y-axis from z=2 to z=0: this ensures that the area under the curves accurately tracks the number of galaxies present. The distributions 
are shown at z = 3.06, 2.07, 0.99, and 0.00. The right panel shows the evolution of the median mass of the same descendants. The median 
of the entire sample is shown in addition to as well as the starforming and quiescent subsamples. It is clear that after z = 3, the median 
mass of the starforming and quiescent galaxies grow at similar rates. 
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Fig. 10. — Left panel; the descendants from the n = 2 X 10~* Mpc~'^ number density selection are shown. Starforming galaxies are blue 
while quiescent galaxies are red The line is the running median for each sample. The points have a star formation floor at log (SFR) = 
-3.7 for illustrative purposes. Right panel: same as the left panel, except now merger growth averaged over time has been added to star 
formation for both galaxy samples. At the highest masses, the mass growth of quiescent galaxies matches that of star forming galaxies. At 
lower stellar masses, the mass growth rate diverges, likely because of the predominance of quenched satellites in this mass range. 



Since velocity dispersion is not reported in the Guo 
et al. (2011) galaxy catalogs, it must be inferred from 
stellar mass and size. We use the following formula to 
calculate stellar velocity dispersion: 



5R 



(1) 



This simple estimate is derived from the scalar virial 
theorem (see, e.g., Cappellari et al. (2006)). is 
the total stellar mass, G is Newton's gravitational con- 
stant, and R is the half-mass radius. The Guo et al. 
(2011) galaxy catalogs provide the exponential scale ra- 
dius of the disk, which is converted into the half-mass ra- 
dius via i?haif-mass = l-678i?exp- The three-dimensional 



half-mass radius of the bulge is converted into a two- 
dimensional projected half-mass radius following work in 
Ciotti (1991) by assuming a,n n — 4 Sersic profile (Sersic 
1968). These two radii are combined into a projected 
half-mass radius of the entire galaxy: 



Rtnt — 



^-^bulge-^hulge 



■disk 



Mhulge + Mdisk 



(2) 



Assuming perfect photometry and a spatially homoge- 
nous stellar population, this would be the observed ef- 
fective radius. 

The upper left panel of Figure [TT] shows the dispersion 
abundance functions as calculated for galaxies from the 
Guo et al. (2011) SAM. They are remarkably similar 
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Fig. 11. — Top left panel: the inferred velocity dispersion function for four redshift slices from the Guo et al. (2011) mock galaxy catalogs. 
The inferred velocity dispersion function is remarkably stable over the whole redshift range. Top right panel: The relationship between 
inferred velocity dispersion and stellar mass at 2 = from the Guo et al. (2011) mock galaxy catalogs. There is no simple one-to-one 
correlation that would make dispersion selection a trivial version of stellar mass selection. Lower left panel: Tracing the evolution of the 
descendants of galaxies from the initial selection at z=3.06. The panels show the distribution of evolution in inferred velocity dispersion 
for individual galaxies. The average is marked with a circle, while the median is marked with a cross. The difference between the 75th 
quartile and the 25th quartile is shown in the upper left of each panel. The n = 5 X 10~^ Mpc~^ descendants are multiplied by a factor of 
two for comparison, to compensate for the smaller initial bin size. Note the stability of inferred velocity dispersion for individual galaxies 
as well as the whole distribution. Lower right panel: The median stellar mass growth as a function of number density and redshift, with 
no adjustments applied. Dispersion selection traces the median stellar mass with remarkable accuracy: compare to Figure |4] 
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across all redshifts. While a stable velocity dispersion 
abundance function implies a stable velocity dispersion 
for individual galaxies, it is not conclusive. The rate 
of change of inferred velocity dispersion for in dividual 
galaxies is shown in the lower left of Figure [TT] as well. 
Comparison to the same plot for stellar mass (Figure Isl 
shows that dispersion is considerably more stable: for 
example, most galaxies change less than 0.06 dex in dis- 
persion from 2; — 1 to z = 0, whereas the spread in 
stellar mass growth in the sam e tim e period is 0.24 dex. 
The upper right panel of Figure [TT| shows the relation be- 
tween stellar mass and velocity dispersion for all galaxies 
in the Millenium Simulation, establishing that the de- 
scribed dispersion selection is unique and not a simple 
function of stellar mass. 

Finally, in the lower right panel, the mass evolution 
of galaxies at a constant cumulative inferred dispersion 
is compared to the mass evolution of descendants of the 
original dispersion selection at z = 3. This is analo- 
gous to Figure |4j except now we are selecting galaxies 
at each snapshot in redshift by their velocity dispersion 
and comparing their stellar masses. Constant cumula- 
tive dispersion selection predicts the median stellar mass 
to within 0.05 dex without correction. The total mass 
evolution of the descendants is 0.6 dex, similar to that 
of the stellar mass selection over the same redshift range. 
This implies dispersion selection traces the stellar mass 
evolution to a much higher accuracy. Also notable is the 
lack of a systematic difference between the inferred and 
true descendants at the higher number densities. There 
is a hint of a systematic discrepancy at the lowest number 
density, but it is small. 

Inferred velocity dispersion is determined from the 
quantity M*/i?. Thus, the applicability of the results 
shown here to the real Universe is dependent on the ac- 
curacy of the evolution of M^/R in the SAM. While the 
stellar mass evolution has been analyzed extensively in 
the previous section, it has been shown that the bulge 
and disk sizes calculated in the Guo et al. (2011) galaxy 
catalogs do not match obser vations. This will be dis- 
cussed further in Section [631 Iii spite of this, the stabil- 
ity of stellar velocity dispersion with respect to merging, 
quenching, and satellites/centrals indicates this may be 
the most promising approach to number density selec- 
tion. 

6. DISCUSSION: CAN NUMBER DENSITY SELECTION BE 
TESTED IN A SAM? 

We have shown that tracking galaxies at a constant cu- 
mulative number density in stellar mass reproduces the 
descendant mass evolution to within 40% over a wide 
range in number densities. Correcting for the scatter 
in growth rates and merging galaxies reduces this er- 
ror to 12%, with the remaining error presumably coming 
from inaccuracies in the correction method and system- 
atics introduced by satellite quenching. The function- 
ality of number density selection and its dependence on 
the merger rate, the quenching mechanism(s), and the 
scatter in stellar mass growth rates are results that are 
independent of the SAM framework. Since the details of 
semi-analytical models and N-body simulations do not 
perfectly mirror the real Universe, however, any quan- 
titative results derived from them must be subject to 
scrutiny before they are applied to observational studies. 



6.1. The Millenium Cosmology 

One of the known issues with the Milleni um Simulation 
is the cosmology: as discussed in Section [Xl] the Mille- 
nium (Tg is nearly 4cr higher than that from the WMAP7 
results (Komatsu 2010). This results in a higher clus- 
tering amplitude for dark matter than is observed in the 
real Universe. The projected autocorrelation function of 
stellar mass for the Guo et al. (2011) galaxy catalogs 
overestimates the clustering amplitude by 10-20 percent 
on large scales and a factor of two on small scales as 
compared to SDSS results from Li & White (2009). 

An inflated erg parameter will affect the stellar mass 
growth rate in two ways. First, it will increase the merger 
rate. Second, it will increase the fraction of quenched 
satellites at a fixed stellar mass. Both effects mean that 
number density selection will more accurately trace the 
median mass evolution of galaxies in the real Universe 
than in the Guo et al. (2011) galaxy catalogs. This is 
explained below for each effect. 

Since mergers are necessarily stochastic events, in- 
creasing the merger rate will increase both the scatter 
in stellar mass growth rates and the fraction of galaxies 
that have " disappeared" in a merger with a larger galaxy. 
If two of the three sources of error are then lower in 
the real Universe, this will only strengthen the efficiency 
of number density selection. Thus, it is expected that 
the discrepancy between the inferred and actual stellar 
mass evolution will be less than 0.15 dex when applied 
to the real Universe. It should be noted that decreasing 
the merger rate will also decrease the mass growth rate 
due to mergers: the remarkable congruence in Figure 10 
between the mass growth of quiescent and starforming 
galaxies may not hold if the merger rate is lowered. 

Satellite quenching is a problem for number density 
selection, as it removes galaxies from the starforming se- 
quence in a stochastic way. Quenched satellites will stop 
growing in stellar mass, whereas a field galaxy of simi- 
lar stellar mass will continue to grow: this scrambles the 
rank order of galaxies. Thus, an increase of the frac- 
tion of quenched satellites will induce more confusion in 
the rank ordering of galaxies. If the real Universe has 
fewer quenched galaxies than the SAM dictates, then it 
is reasonable to expect that number density selection will 
trace the median stellar mass with greater accuracy than 
in the Guo et al. (2011) SAM. 

It is also noteworthy that the mass growth rate will 
affect the redshift evolution of the stellar mass function. 
Since the Millenium cosmology results in many galax- 
ies quenching too "early", the semi-analytical model is 
forced to build up the low-mass end of the galaxy mass 
function at early times. This results in an excess of low- 
mass galaxies at high redshift compared to observations- 
see Guo et al. (2011) for in-depth discussion. This dis- 
crepancy between the semi-analytical mass function and 
the observed mass function is a result of a prescription 
for [M^ ) that differs from the real Universe. At least 
some of the issues causing this discrepancy have been 
addressed above; however, until semi-analytical models 
can accurately produce the redshift evolution of the mass 
function, their quantitative predictions for galaxy stellar 
mass evolution cannot be fully trusted. 

6.2. Mass Growth of Quiescent Galaxies 
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In order for the number density selection scheme to 
be successful, quiescent galaxies must experience signifi- 
cant mass growth through galaxy-galaxy mer gers. This 
is true in the SAM, as evidenced by Figure [TOl Evidence 
for the importance of mergers in growing stellar mass 
comes from studies of cosmological and N-body simu- 
lations (Naab, Johansson, & Ostriker 2009; Oser et al. 
2010; Feldmann et al. 2010; Lackner et al. 2012). Ad- 
ditionally, the tidal tails and streams seen in deep imag- 
ing of local massive quiescent galaxies support the pic- 
ture of mass growth via mergers (van Dokkum 2005; Tal 
et al. 2009; Janowiecki et al. 2010). Thus, it is a reason- 
able expectation that quiescent galaxies still experience 
a significant amount of stellar mass growth via merg- 
ers. Whether it is sufficient to match the mass growth 
of starforming galaxies- as is true in the SAM- requires 
a careful observational study that is beyond the scope of 
this paper. 

6.3. Size Growth and Inferred Dispersion 

The efficiency of ordering galaxies by their velocity dis- 
persion, as opposed to their stellar mass (see Section 5), 
is dependent upon conservation of rank order in Af^/i?. 
This means that accurately modeling the size evolution 
of galaxies within the SAM is important. In the semi- 
analytical model, disk sizes are set by the angular mo- 
mentum of their stellar component, which is in turn set 
by the angular momentum of the gaseous component at 
the moment it is transformed into stars. The specific an- 
gular momentum of gas accreted from the IGM is set to 
be equal to the specific angular momentum of the halo; 
gas acquired during mergers has its own angular mo- 
mentum which is added to the galactic gas. Bulges grow 
during major mergers, minor mergers, and due to disk 
instabilities. Their size growth is determined by energy 
conservation and the virial theorem. 

Guo et al. (2011) notes there are systematic differences 
between the resulting distribution of bulge and disk sizes 
as compared to observations. Both the median bulge and 
the median disk sizes are too high at low stellar masses, 
and too low at high stellar masses. Additionally, the 
simulation scatter is larger than the observed scatter. 
In particular, the small observed scatter in stellar bulge 
size with stellar mass is difficult to understand if bulges 
are built primarily via mergers (Nair, van den Bergh, & 
Abraham 2010), as is true for high and low mass galaxies 
(although not intermediate. Milky Way mass galaxies) in 
the Guo et al. (2011) model. This casts doubt on the ac- 
curacy of stellar mass predictions derived from ordering 
galaxies by their velocity dispersion. 

However, the current model for size evolution incor- 
porates most known important physical processes and 
should not be dismissed out-of-hand. For example, the 
velocity dispersion function is believed to be remarkably 
stable (Bezanson et al. 2011), and this stability is re- 
produced in the model. Additionally, as reflected in 
the merger model in Guo et al. (2011), some physical 
processes are expected to couple mass and size growth. 
In this scenario, velocity dispersion may remain a rela- 
tively stable quantity. Exploring whether constant num- 
ber density in velocity dispersion, in stellar mass, or sim- 
ple constant velocity dispersion is the best approach to 
tracking galaxies will require a more sophisticated and 
robust semi-analytical recipe for velocity dispersion, pos- 



sibly attainable by comparison to hydrodynamical simu- 
lations. 

6.4. What Results are Independent of SAM Recipes? 

There are certain features that must exist in a semi- 
analytical model in order for it to be a reasonable testbed 
for number density selection. The star formation rate 
must growth with the stellar mass of a galaxy, and the 
merger rate must increase monotonically as a function 
of stellar mass. A quenching model that turns off both 
massive centrals and satellite galaxies must exist, in par- 
allel with the real Universe. Spontaneous galaxy creation 
must be a negligible effect at higher stellar mass. These 
features are all evident in the Guo et al. (2011) SAM. 
This allows us to extract conclusions from the behavior 
of galaxies in the SAM which may then be applied to 
observational studies, independent of details in the semi- 
analytical and numerical modeling. 

The simple techniques used to adjust the inferred mass 
growth rely only on a few observables: the shape of the 
stellar mass function, the merger rate as a function of 
stellar mass, and quantification of the starforming se- 
quence. All of these may be estimated from observations, 
although some for example, the merger rate are poorly 
constrained, particularly at high redshift (see, e.g., New- 
man et al. (2012)). Thus, even if the specifics of the semi- 
analytical model prove to be incorrect, number density 
selection will still operate within the same framework as 
described here, making the method can easily generaliz- 
able to the real Universe. 

The numerical results presented here remain sensitive 
to the details of the SAM recipes. The slope, scatter, 
and exact details of the merger rate and the starforming 
sequence affect the mass evolution as a function of num- 
ber density. The specifics of the quenching recipe may 
also prove critical to the efficacy of tracking galaxies at a 
constant number density. With reasonable prescriptions 
for these effects, though, the resulting inferred and de- 
scendant mass evolution should be good estimates of the 
truth. 

6.5. Comparison to Tracking Galaxies with a Stellar 

Mass-Halo Mass Relationship 

The evolution of halo mass is known to high accuracy 
from numerical simulations. By adopting a redshift- 
dependent stellar mass-halo mass relationship, one can 
use the evolution of the halo mass function to infer the 
evolution of stellar mass (e.g. (Conroy & Wechsler 2009; 
Behroozi, Wechsler, & Conroy 2012)). The strength of 
this approach relative to number density selection is that 
the growth and merger rates of dark matter halos are 
known to greater accuracy than the growth and merger 
rates of galaxies. However, an additional source of scatter 
is introduced when the stellar mass-halo mass relation- 
ship is applied. This scatter is on the order of 0.16-0.2 
dex at z = 0, with the scatter presumably increasing 
with redshift (More et al. 2009; Reddick et al. 2012). 
A comparable analysis to that presented in this paper 
would be needed is necessary to determine the effects of 
the scatter in the stellar mass-halo mass relationship on 
the inferred stellar mass evolution; unfortunately, studies 
parametrizing M*(M/i,z) separately for individual halos 
have not yet been done (Behroozi et al. 2012). Pending 
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further stiidy, then, tracking galaxies with information 
from dark matter simulations may prove to be a viable 
alternative to number density selection in the future. 

7. CONCLUSIONS 

We have investigated the efficacy of tracking galaxies 
at a constant cumulative number density. We discuss the 
assumptions inherent in the use of number density selec- 
tion and which physical processes might violate them. 
We then use the Guo et al. (2011) semi-analytical model 
based on the Millenium Simulation to examine number 
density selection in a realistic simulation setting. The 
technique is demonstrated to be a robust method to track 
the evolution of an individual galaxy population. Apply- 
ing the technique to galaxies evolving from ^ = 3 to 
2; = 0, or eleven billion years in time, reproduces the 
factor of four evolution in stellar mass to within 40% 
within the model. Effects that confuse stellar mass rank 
order have been identified: scatter in stellar mass growth 
rates, merging galaxies, and quenching. Corrections de- 
veloped for the first two effects brings the discrepancies 
down to approximately 12%. However, a difference in 
the mass growth rate of quenched galaxies relative to 
starforming galaxies may produce systematic errors in 
the assumption of constant rank order. This is sensitive 
to both the recipes for quenching and the merger rate 
and thus merits investigation in the real Universe rather 
than within a semi-analytical model. We also examine 
ordering galaxies by their inferred velocity dispersion as 
opposed to their stellar mass. This method reproduces 
the mass evolution of galaxies from z = 3 to 2 = to 
within 12% with no corrections whatsoever, although it 
is sensitive to the rules governing size evolution in the 
semi-analytical model. 

The ultimate goal of this study is to apply the number 
density technique to observational studies of the real Uni- 
verse. However, since it is impossible to observe a galaxy 
evolve through time, we must calibrate this technique in 
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semi-analytical models. The most obvious concern in this 
study is that semi-analytical models may not accurately 
reflect the Universe. Some of the most important discrep- 
ancies include the redshift evolution of the mass func- 
tion and enhanced clustering strength compared to the 
real Universe. We address these concerns in a number of 
ways. First; we emphasize the SAM is used primarily as 
an investigative tool to see how number density selection 
works and what will disrupt it: the major mechanisms 
(mergers, the star forming sequence, quenching) affecting 
this are reproduced in the SAM, regardless of the details 
of their implementation. Second, the methods developed 
herein to correct for scatter in growth rates and mergers 
arc based on simple observable quantities like the stel- 
lar mass function and the merger rate, and may safely 
be generalized to the real Universe. Finally, it is shown 
that the known discrepancies between the semi-analytical 
model and the Universe act to decrease the errors in the 
observational inferred mass evolution relative to those in 
the SAM. Thus, the errors quoted here should be taken 
as upper limits to errors that would arise in an analogous 
study of the real Universe. 

Linking descendant and progenitor galaxies at a con- 
stant number density is thus shown to be a considerable 
improvement over previous techniques, and opens the ex- 
citing potential for using observations to track changes in 
individual galaxy populations over long periods of cosmic 
time. 
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